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Abstract. Non-unitary operations generated by an effective non-Hermitian 
Hamiltonian can be used to create quantum state manipulations which are 
impossible in Hermitian quantum mechanics. These operations include state 
preparation (or cooling) and non-orthogonal state discrimination. In this work 
we put a lower bound on the resources needed for the construction of some given 
non-unitary evolution. Passive systems are studied in detail and a general feature 
of such a system is derived. After interpreting our results using the singular value 
decomposition, several examples are studied analytically. In particular, we put 
a lower bound on the resources needed for non-Hermitian state preparation and 
non-orthogonal state discrimination. 
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1. Introduction 



Unitary quantum evolution generated by Hermitian Hamiltonian conserves the angle 
between states. For example, two states which are orthogonal at one instant will 
remain orthogonal throughout the evolution. In effective models of some sub-systems 
the Hamiltonian may be non-Hermitian. This may be either due to leakage of particles 
outside the region of interest, (as happens, for example, in Stark ionization), or due 
to absorption or gain, as often occurs in optics. Another source of non-unitarity can 
appear when some measurements are performed on entangled systems (e.g [I]). For 
additional information see Effective non-Hermitian Hamiltonians lead to non- 
unitary evolution operators. It was suggested in [3] that non-Hermitian Hamiltonian 
can be used for perfect discrimination of two non-orthogonal states. In principle this 
scheme can work even if the states are almost perfectly parallel to each other. In [3] 
the question of the needed resources for this task was not addressed. In this work, we 
have fixed a lower bound on the resources needed for a non-unitary task in order to 
address this question. 

We start in Sec. [Uby introducing the "norm action" as a plausible measure for the 
resources needed to generate some evolution. In Sec. |3]we derive our basic result: a 
lower bound on the norm action that depends only on the evolution operator at the end 
of the process. That is, a simple bound can be set on the desired evolution resources 
without knowing anything on the Hamiltonian. Passive systems are discussed in sec. 
[¥] where our basic result is generalized so that it will be effective for passive systems 
as well. In Sec. [5] we use the singular value decomposition to explain in what sense 
our bound quantifies solely the "purely non- unitary" part of the evolution. Next, in 
section |6] we analytically study several examples in order to illustrate various points 
concerning our main result. Section [7] contains concluding remarks and the three 
appendices contain some clarification. 



2. Quantifying evolution resources 

Let U{t) G C^^^ be some some non- unitary evolution operator in an A^-level quantum 
system. The action of U{t) on a state \tp{t = 0)) is: 

\m) = uit) \m) ■ (1) 

In general, the determinant of U{t) is not fixed in time, and it satisfies the Liouville 
formula ^: 

det[i7(t)] =det[L/(0)]e-*''"*''[^^*')l'^*', (2) 

where H is the Hamiltonian of the system. The trace of H only contributes a complex 
multiplier to the evolution operator. Obviously the trace affects the required resources, 
but these resources are rather trivial as they amount to adding a common phase 
or attenuation to all states. In the first stage we are only interested in resources 
associated with actual state transformation. Therefore, to avoid the trivial changes in 
U generated by the trace of H we define the traceless Hamiltonian: 

n = H-^i. (3) 

The evolution operator U generated by this Hamiltonian satisfies: 

det[U{t)] = 1, (4) 
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and it is related to U through: 

U = U/det[U]^/^ (5) 

U is the determinant normahzed version of U. The geometric meaning of condition (|4]) 
is given in appendix I. Despite note that U is in general not a unitary operator: 

U'^U ^ I, (6) 

where f stands for the complex transposition operation. Assume there is some non- 
Hermitian Hamiltonian H{t) that after some time T generates a desired evolution 
U(T). Our goal will be to use U{T) to set limits on the smallness of H. However, 
a simple scaling argument shows that for any final U, the generating Hamiltonian 
elements can be made arbitrarily small by making the evolution slower. To remove 
this ambiguity we will examine the time integral over the Hamiltonian "size". The 
scalar quantity we shall use to quantify the Hamiltonian size at each instant is the 
matrix norm of the Hamiltonian matrix. In summary, our goal is to set a lower bound 
on the action-like quantity: 

"^Wnmdt (7) 







We shall refer to this quantity as the "norm action". In this work ||-|| refers to any 
matrix norm that satisfies the sub-multiplicative property: 

\\AB\\<\\A\\\\B\\, (8) 

In particular, we shall often use the spectral norm [S] of the Hamiltonian. A 
previous study [6 showed that the spectral norm of Ji constitutes an upper bound 
on the evolution speecQ of a quantum state driven by NH or Hermitian Hamiltonian. 
This physical meaning of the spectral norm makes it our default choice of norm. 
Nonetheless, some of the results presented here equally apply to any sub-multiplicative 
norm just as well. Therefore the spectral norm will be used only where it can provide 
some additional insights unique to this norm. The (sub-multiplicative) spectral norm 
of the Hamiltonian is the largest singular value of H [7] and it is equal to: 



\\H\\^ = ^max[eig{i/ti/}], (9) 

where eig{} is just the list of eigenvalues of the enclosed operator. Before we proceed, it 
is worth mentioning two properties of the spectral norm which provide more explicit 
information on the Hamiltonian matrix elements. Let Hij be the elements of an 
operator H £ C^^^ . The spectral norm satisfies: 



N 



max|F,j| < \\H{t)\\^ < Nmax\H,j\ (11) 
That is, a lower bound on ([7]) would also constitute a lower bound on 



/o \/Y.t,j=i\Hij{t)\ and on max\H,^j{t)\ dt. 

X Loosely speaking the angular speed in "Bloch Sphere" 
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3. A lower bound on the resources needed for non-unitary evolution 

In this section we derive a lower bound on the "norm action" ([7]). The derivation is 
carried for the more general H and U and only later replace them by H and U. The 
evolution operator U{t) satisfies: 



i — U = HU 
dt 



Taking the norm and using sub-multiplicativity 



dt 



or: 



li^li>llt/|| 



dt 



(12) 



(13) 



(14) 



Since any norm satisfies the "inverse triangle inequality": HA — _B|| > — ||-B|||, we 
get: 



dt 



> 



\U{t + dt) - U{t)\\ 
Jt 

\U{t + dt)\\^\\U{t)\\ 



lim 



Therefore: 



\H\\>\\U\\ 



dt 



dt 



\u{t)\\ 



dt 



\u{t)\\ 



>\\ur'j,\\um 



Integrating we get our basic resources bound: 



\H{t)\\ dt > \n\\U{T)\\ 



(15) 



(16) 



(17) 



An alternative derivation is given in appendix II. 

Although this equation is correct, it is sensitive to tr(iJ) and det[[/] which are 
related to the somewhat irrelevant information discussed earlier. This can be easily 
resolved by noticing that Ti and U also satisfy i^U = HU, so pT|) can be used to 
obtain: 



\nit)\\dt>\n\\U{T)\\ 



(18) 



One may be concerned that perhaps ||Z^(T)|| < 1 and therefore the RHS of (|18p is 
negative which makes Eq. ((T8|) trivial and useless. Nonetheless, for the spectral norm, 
one can show that: 

||Z^(r) 11^ = 1 is unitary (19) 

Thus, for non-unitary operators the RHS of (|18|) is larger than zero and our lower 
bound is non-trivial. A simple derivation (|19p is given in appendix III. Using the 
following property of the spectral norm: 

||A||^ > max|eig{A}| (20) 
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For any operator A e C^^^, it is possible to obtain a weaker lower bound in terms 
of the eigenvalues of evolution operator: 



In our third example we will examine a non-diagonalizable case (NH degeneracy) 
where all the eigenvalues are equal to one so that the RHS of (|2T|) is equal to zero. 
Nonetheless, the operator is still non-unitary and the stronger resources lower bound 
(|18p is greater than zero. 

The fact that the RHS of ([T5)) is equal to zero for unitary evolution operators (due 
to and ([TO]) ) gives some indication that our spectral norm bound quantifies the 
resources needed for non-unitary evolution but not the resources for unitary evolution. 
This point will be clarified in Sec. [SJ However, before doing that we wish to extend 
our result to better fit the common case where the trace of the Hamiltonian is non-zero 
(e.g. a dissipative system). 

4. Improved inequality and the marginally passive case 

While losses are very common in physical systems, gain requires some special devices, 
special materials, power sources, etc. Therefore it is often easier to implement non- 
unitary evolution in passive systems (zero gain) at the expense of smaller signals at 
the output of the system (e.g. see the discussion in [9]). As will be explained next, the 
basic result (fT7|) . does not provide any useful information for passive systems as the 
inequality is trivially satisfied. The goal of this section is to generalize the previous 
results and remedy this problem. 

Mathematically, a passive system cannot increase the norm of an input state. At 
most it can preserve it. In principle there is no point in attenuating the signal more 
than needed to reach passiveness. We shall use the term "marginally passive evolution" 
to describe a passive system that has at least one state whose norm is unchanged by 
[/, while the norm of other states is either left unchanged or reduced. By definition the 
maximal norm of an output vector is given by the spectral norm of the transformation 
U. The final amplitude of a state {ipj) generated by U from a normalized initial state 
lipi) is simply: ^ (V'/l'0/) — V {''Pi \U^U\ ipi) . In fact, an alternative definition of the 
spectral norm is given by: 



Therefore, given a general non-unitary U, the marginally passive (MP) evolution 
operator is given by: 



Since ||C/mp|Is = 1, and more generally for passive evolution || C/pH^ < 1, it becomes 
obvious that (|17p is trivially satisfied for any passive evolution (LHS is positive while 
the RHS is negative or zero) . In what follows we improve (|17p to account for the 
passive system as well. 

Inspection of (fT7|) reveals that there is an unnatural asymmetry in the RHS. The 
"magnitude" of the non-unitary operation is determined by the maximal amplification 
(maximal singular value) while the maximal attenuation is completely ignored. 




(21) 




(22) 
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Although in practice amplification is much more difficult to set up in comparison to 
attenuation, from a mathematical point of view, they should be treated on an equal 
footing. This asymmetry manifests itself clearly in the MP case. The amplification 
is zero and all the non-unitary information is in the degree of attenuation of state. 
To extract the attenuation information we will consider the inverse operation U~^. 
The logic behind doing so, is that amplitude attenuation when going forward in time 
is equivalent to amplification when going backwards in time. We start by taking the 
time derivative oiU^^U = I form which we get: 

dtU = -UidtU-')U (24) 

upon using it in the Schrodinger Eq. ([T2|) we get: 

- idtU-^ = U-^H (25) 

Taking the norm on both sides and using sub-multiplicativity as we did before, we 
get: 

||atC/-^|| < ||[/"^|| (26) 
Repeating the same procedure used to derive (jl7p we obtain: 

rT 

\\H\\dt>\n\\U-'^\\ (27) 







Combining it with (jl7p we get the main result of this paper: 



T 



\H\\ dt > ln[max{||J7|| , (28) 



Since ||[/||^ ||C/ -'-||^ > 1 , for passive systems ||C^p"'"||^ > l/l|t^p|ls ^ 1- If the 
attenuation in the system provides a nontrivial result via the term. Equation 
(|28p resolves the aforementioned asymmetry between gain and loss. Later on it will 
be shown more explicitly for the spectral norm. 

There is another way to combine the ((T7|) and ((27)) which is less efficient than 
(|28l) (i.e. a smaller lower bound), but it offers some insights into the interplay of gain 
and loss and a more analytical structure. By adding ([T7|) and (f27ll and dividing by 
two we get: 

''^\\H\\dt>\nVmW^\ (29) 







The norm of U in p7p is replaced by the geometric mean of the norms of U and U ^ . 
Moreover this expression is invariant under changes of the determinant of U (or trH) 
since detU = 1/detL/"^ 

The results (pS)) and apply to any sub-multiplicative norm. Now we wish to 
focus on the spectral norm. It is easy to show that the singular values, s, of U, are just 
the inverse of the singular values of U^^ (provided that U is invertible). Therefore: 

/ \\H\\^dt > max{ s„a:., s;;^„} > In J s^axs:;;,]^ (30) 

Jo 

where Smax and Smin are the largest and smallest singular values of U. Here the 
symmetry between gain and loss appears more explicitly. 
In the marginally passive case ||[/a/p||^ = 1 so: 

\\HMPtdt>\n\\Uj^/pt (31) 
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Finally, to conclude our discussion on the passive case we would like to point 
out a general feature of passive systems. It concerns the tradeoff between "non- 
Unitary performance" and "detectability". As shown in appendix I, |det[/| can 
be interpreted as the geometric average of the state amplifications generated by U. 
Starting from Ump and using ^ we get: 



WUL. = 



l/N 



UMp/detpMP^ 

||C/A/p||^/|det[/MP 
l/|detC/Mp|'/^ 



l/N 



(32) 



Combining this with the expression for the geometric mean amplitude, ((a)), from 
appendix I we get: 

1 



{{aMp)) = det[UMp]'/'' 



ll^ll. 



< 1 



This clearly shows the performance vs. detectability tradeoff in the marginally passive 
case. The "stronger" the non-unitary operation, is the bigger \\U\\ is. Yet, as \\U\\ goes 
up, the geometric mean amplification goes down. In other words, if the system has no 
gain then stronger non-unitary operation will lead to a weaker signals. 



5. "Purely non-unitary" evolution and unitarily invariant norms 

In [5] a somewhat similar formula to psp was written for unitary evolution. It follows 
from [8] that in the unitary case our In ||Z//(r)|| is replaced by ||lnW(T)||. Unfortunately 
the appealing derivation in [8] is suited only for Hermitian Hamiltonians and the 
analysis cannot be easily extended to NH Hamiltonians and non-unitary evolution. 
Our derivation overcomes this difficulty (also up until now we did not require that the 
norm be unitarily invariant), at the expense of the capability to capture the resources 
associated with unitary evolution. This can actually be an advantage if indeed our 
result always captures only the resources needed for the "purely non-unitary" part of 
the evolution. In this section we clarify in what sense and under what conditions our 
result (|18p quantifies only the purely non- unitary part of the evolution operator. 

For simplicity let us assume that ZY(T) is diagonalizable and that {fii} are the 
eigenvalues of the evolution operator. If all eigenvalues satisfy \fii\ = 1, as in the 
unitary case, the information on the evolution is encapsulated in Arg[fii]. However 
if the eigenvalues are real and positive (and not all of them equal to one) , it is clear 
the evolution is "purely non-unitary". The eigenstates are just attenuated or amplified 
without accumulating any phase as in the unitary case. With this in mind, we can 
make use of the singular value decomposition (SVD). According to the SVD scheme, 
any operator U{T) e £,nxn (■gygj^ jf jg j^q^^ diagonalizable) can be decomposed into 

HI: 

U{T) = QWP^ (33) 

where Q and P are some unitary matrices!^ and is a diagonal non-negative matrix 
whose entries are the singular values mentioned earlier. In the unitary case W is the 
identity matrix. In the non-unitary case W corresponds to a purely non-unitary 
evolution operator. Equation (|2ip shows that any evolution can be written as a 



§ The exact forms of P and Q are not needed for our purpose. 
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product of unitary evolution followed by a purely non-unitary evolution and then 
by another unitary evolution. Until now we did not require that the norm be unitarily 
invariant (UI). However, unitary invariance facilitates a very clear interpretation of 
the results. A norm is UI if it satisfies: 

\\A\\^, = \\CAD\\^, (34) 

for any A G C^^^ and any unitary C and D. The spectral norm, for example, is 
unitarily invariant norm. Using this property and Eq. p3p we get : 

\n\\U{T)\\^j^\n\\W\\^j (35) 

That is, the unitary part of evolution is removed and the lower bound we found 
addresses only the purely non-unitary part of the evolution. 



6. Examples 

In this section we demonstrate a few points about the resources inequality through 
explicit examples. 



6.1. Example 1 - A simple decay 

Consider a two-level time- independent Hamiltonian with one decaying level and one 
stable level: 



^ = ( e ) 

The traceless Hamiltonian and the corresponding evolution operators are 



U=[^ „ir* ) (38) 

iir 

-¥ 



'I ) (39) 



The norm action is: 

nT 







H{t)\\^dt=\vT (40) 



The RHS of the resources inequality yield: 

\n\\U\\^ = \TT (41) 

This elementary example demonstrates two important points. The first, is that the 
resources inequality (|18p can be an exact equality in some cases. The resources 
inequality, then, is not an exaggerated overestimate. The second point is the 
importance of the trace removal and the U normalization. For F > it is easy 
to see that \\U\\^ = 1 so ln||J7||^ ~ which make the resources inequality useless 
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(/o ll-'^o(i)|ls > 0). Alternatively. F > is exactly the marginally passive case so 

(42) 
(43) 



instead of normalizing U we can use (|28p and get: 

[ WHom^dt^TT 
Jo 

Which is again an equality rather than an inequality. 
6.2. Example 2 - A cooling evolution 

Consider a quantum two-level system that has two solutions: 



IV-i) 



1 



sin 



cos 



m 

4 



(44) 
(45) 



where 9{0) ~ 0. For simplicity we consider 9'{t) > so that 9{t) is a monotonously 
increasing function. While the two solutions are initially orthogonal to each other, 
as — ^ TT the two states become parallel to each other. This transformation changes 
a mixed state density matrix into a coherent one. Alternatively, this scheme can be 
used for non-orthogonal state discrimination. If one starts at 6'(0) 7^ and decreases 
the angle 6, two initially non-orthogonal states become orthogonal when 6* = 0. Once 
the states are orthogonal they can be resolved with 100% certainty. This kind of 
non-Hermitian state discrimination was first proposed in [3J. However, in [3] the 
needed resources were not considered. Moreover, in our example the Hamiltonian is 
time-dependent. The two solutions p4)) . (|45|) can be used to construct the evolution 
operator U from which U and H can be explicitly calculated: 



n 



U = 




\\u\\ = 



' 1 + sin ■ 



cos ; 





1 


1 

2 


Tan 


e(t)l 




2 




e{t) 


) 


2 





2 



Now let us evaluate the norm action. The Hamiltonian norm is: 

4 cos[2] 
assuming monotonicity {9' > 0): 



sm - 



(46) 
(47) 

(48) 
(49) 

(50) 



Resources needed for non-unitary quantum operations 



10 



After using few trigonometric identities we find that: 

ln||W|L = xln — = -\n — (51) 

2 cos I ^ COS I — sm I 

Equations ((50l) and ((5T|) show that the resources inequahty ((T8|) holds {9 > 0). In this 
example the resources inequahty was demonstrated for time-dependent HamiUonians. 

The resources divergence at 6* tt is a typical feature of non-unitary coohng and 
it can be understood without a detailed analysis. Since at the end of the evolution 
the two states become almost parallel to each other, the column vectors of U become 
linearly dependent. As a result det[J7] ^0. In order to to keep det[Z//] = 1 (see ) the 
element of U must be very large. The LHS of Eq. (|TT|) shows that if the elements of 
an operator are very large, the spectral norm of the operator is also very large. 

The divergence at 6* vr seems to indicate that perfect cooling is not possible. 
Nonetheless, the divergence is only logarithmic, so in practice the norm action does 
not obtain large values in this case. For example, if 6* = tt — lO"'^^, the minimal 
required resources is only about 15.3 in units of h. 

Nonetheless, the divergence is only logarithmic so in practice the norm action 
does not obtain large values in this case. For example, ii 9 = w — 10~^^ the minimal 
required resources is only about 15.3 in units of h. 

The comparison of NH cooling or state discrimination to other methods 
(relaxation via the Lindblad equation, POVM) is outside the scope of this work. Our 
goal is only to investigate fundamental limitation on the construction of non-unitary 
operations generated by the Schrodinger equation with NH Hamiltonians. 



6.3. Example 3 - Non-Hermitian degeneracy 

In this example we wish to examine the resources in the most extreme non-Hermitian 
scenario - the exceptional point Hamiltonian. Exceptional point Hamiltonians are 
Hamiltonians that cannot be diagonalized. At most, similarity transformations can 
bring these Hamiltonians to the standard canonical Jordan Block form. See Ch. 9 of 
(2] and references therein for more information on exceptional points. Consider the 
following exceptional point Hamiltonian and its corresponding evolution operator: 



(52) 

The norm action is: 

f \\H{t)\\^dt^\Eo\T (54) 
Jo 

Since is just U with Eq — > — E'o, the bound based on gives the same result 
as the bound based on U . The RHS of the resources inequality is: 



In \\U\l = In ^ 1 + i |i?o| T (|£;o| T+^A+ \Eo\' T^) (55) 
If the resources inequality holds we should expect: 



e^lSoiT >l + ^ \Eo\T{\Eo\T+^^+\Eo\'t^), (56) 
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Next we wish to show analytically that this inequality indeed holds. 



^l + ]^\Eo\T{\E^\T + 2+\Eo\T) 




> 1 + i |i?o| T (|i?o| T + yiTST?^) 



(57) 



This example shows that in some cases the eigenvalues of U are completely inadequate 
to describe the dynamics. In this example eig{Z^} = {1,1} so the RHS of (PT|) is 
equal to zero and therefore the eigenvalues of lA carry no new or useful information 
on the resources integral. This is similar to the finding in [6j, the evolution speed 
near NH degeneracies in the Hamiltonian is determined by the spectral norm of the 
Hamiltonian and not by the eigenvalues of the Hamiltonian. Moreover this example 
shows that time- independent Hamiltonian does not necessarily have a norm action 
which is equal to the bound as was the case in our first example. 

7. Concluding remarks 

It was shown that non-unitary evolution requires some minimal resources. The "norm 
action" of the Hamiltonian has to be larger than the norm of the evolution operator 
or its inverse. For a passive system it was also shown that larger non-unitary changes 
inevitably lead to weaker signals. 

The spectral norm of the Hamiltonian was studied before in the context of evolution 
speed in non- unitary evolution [B]. Moreover it was used to find the most efficient 
Hamiltonians for a given trajectory in the projective Hilbert space. Will these efficient 
Hamiltonian also give the minimal norm action ([7])? Not necessarily. The set of 
constraints in this work is completely different. Here we consider the final state of all 
possible initial states while ^ concerns the resources optimization of a single specific 
state with a given evolution. Furthermore, in this work there is no constraint on the 
evolution between the beginning and the end of the evolution. 

Our bound captures only the "pure non-unitary" part of the evolution. In [8J a norm 
action bound for unitary evolution and Hermitian Hamiltonians was introduced. It 
would be interesting to find a bound that exclusively captures the unitary part of 
the evolution (for non-unitaty evolution) as well as a bound that captures everything 
together and study the relations between the different bounds. 
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Appendix I - The geometric mean amplification 

We define the mean amplification in the following way: let {\'m)}^^^ be an orthogonal 
set under the standard inner product. The amplification of each element generated 
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by an evolution operator U is given by: 



{m\mU\m) 



{m\m) 



The geometric mean is: 



Oavg = V"i---"JV = (li^ eiglU'fU}) = |detC/|'^" 



(58) 



(59) 



To illustrate the reasoning behind the geometric mean consider a two-level system 
with two orthogonal eigenstates. One is amplified by a factor of two and the other by 
a factor of one- half. The algebraic mean is equal to 5/4 while the geometric mean is 
equal to one. The geometric mean, then, matches our expectation that the two and 
one-half should balance each other. If so, the requirement |detZ//| = 1 (or equivalently 
that jj^ tr{H)dt' — 0) used in section (|3]), actually means that the geometric mean 
amplification is equal to one. 

Appendix II 

Given some time dependent Hamiltonian H{t) the evolution operator at some time T 
can be written as: 



U(T) = lim e"*-^'^)''* e~^H(tj)dt ^-iH{u)dt^-iH{to)dt 



(60) 



where tj = j dt. Using the sub-multiplicativity and the basic norm properties we get: 



\\U{T)\\<\im 
Moreover, each term satisfies: 



-iH{T)dt 



-iH{tj)dt 



-iH(to)dt 



(61) 



-iH{tj)dt 



E 



i~zH{t,)dt)' 



k=0 



k\ 



o° Wi-iHUAdty 



k=0 



k\ 



< 



fe=0 



\\dt 



kl 



Therefore: 



\\U(T)\\ < lim ell^^^)"''* g||i/(t,)||dt ^\\H{to}\\dt 
~ dt^o 



Now the RHS is a multiplication of numbers so: 

||[/(r)|| < e'''"'"^"^^""'^^*^^"'^* = e'™"*'^" ^j"^^*^'^"'** 
From which we once again obtain the resources inequality: 



[ > ln||C/(T)|| 

Jo 



(62) 
(63) 

(64) 
(65) 

(66) 
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Appendix III 

In this appendix, for completeness, we give a simple derivation (|19p : 

\\U{T)\\^ = 1 is unitary 

One direction is trivial: if U is unitary so that WU = I, then all the singular 
values are equal to one. In particular the largest singular value (the spectral norm) is 
equal to one. 

For the other direction we use the fact that detU = 1. From this we get that: 

detU^ = {detU^y = 1 (67) 

From which follows: 

Y[sf = detU^'U = detU'' detU = 1 (68) 

where Si > are the singular values of U (the square root of eig{U^U}). If 
||Z//(T)||^ — maxjsi} = 1 then to satisfy it must be that all the singular value are 
equal to one. Thus, the SVD decomposition of the U is: 

U = giypt ^ Qpt (69) 
where Q and P are unitary matrices and therefore: 

U'fU ^ PQ^QP'' = / (70) 
which completes the proof. 
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